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Abstract. I consider production of heavy quarks in pA and AA collisions 
in the framework of the Color Glass Condensate. I discuss the heavy quark 
production in pA collisions in a classical approximation which takes into account 
all multiple rescatterings of a proton in a nucleus. For peripheral collisions, heavy 
quark production cross section can be written in fcy-factorized form. The k^p 
factorization is used to construct a simple model which takes into account both 
classical and quantum effects in pA and A A collisions. I review the main result 
of calculation based on that model: open charm in the central rapidity region at 
RHIC gets suppressed as a function of rapidity. Although the numbers obtained 
for the suppression factor are model dependent, the very fact of suppression is 
the general feature of the Color Glass Condensate at RHIC kinematical region. 
It indicates the onset of quantum evolution effects in a nucleus and is universal 
for all inclusive processes. 

The process of heavy quark production in pA collisions at high energies has three 
separated in time stages in the nucleus rest frame. Emission of a gluon g by a proton's 
valence quark q v takes much longer time T q ^_, qvg than a subsequent emission of a qq 
pair by a gluon t g ^ q q -C i~ qv ^ qvg . In turn, the time of interaction of a q v gqq system 
with a nucleus is of the order of nuclear length Ra and is negligible as compared to 
the evolution time of the proton wave function T qv —>q vg 3> tg^qq 3> Ra- For example, 
the gluon emission takes time 

2q+ 1 



q 2 xMn 



>i? A ~r in t, (1) 



where the proton is moving in + light-cone direction. The same argument applies 
to the successive emission of a quark-antiquark pair by a gluon in a proton's wave 
function. The separation of times is ensured by strong ordering q + 3> k + 3> kx, where 
k + and k^ are the quark energy and transverse momentum respectively. Therefore, 
processes in which gluon or heavy quarks are produced in course of the rescatterings 
in a nucleus are suppressed by powers of energy p + . 

As the result, the differential cross section for the quark production can be written 
in a factorized form as a convolution of the valence quark and gluon wave functions 
with the rescattering factor. The final results takes form [3 El El 

da f j2 u f ,2 f j f d 2 xd 2 y 
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where I assumed for simplicity that the dominant contribution comes from interaction 
of the qq pair with the target, while rescatterings of the gluon and the valence quark 
are neglected. However in general, they must be taken into account as well. In Eq. © 
I used the following notations: 

a s C F ("2 + (1 - a )x ) ■ {ay + (1 - a)x ) 

(ax + (1 - a)a: ) 2 ( a U + (1 - a feo) 2 



7T V k-x \\y-x \ 



x[a 2 + ^o(k-al m )^o(|y ~x \m) J , (4) 

where Q s is the saturation scale, to is a quark mass and a = k + /q + . 

In the nucleus light-cone frame Eq. J5J can be interpreted as a cross section for 
proton scattering off the qq fluctuation of the strong nuclear gluonic field (Color Glass 
Condensate). At the classical level, this field is of the order F ~ Q 2 s /g- It is well known 
in QED that strong electric field E can produce quark-antiquark pairs from vacuum 
if change of its energy eE along distance Ti/mc is larger than toc 2 i. e. eE > to 2 . 
Similarly, strong chromoelectric field can produce pairs from vacuum if Q 2 > to 2 
jS] . Theory of Color Glass Condensate predicts dependence of the saturation scale on 
energy and atomic number: Q s oc A x ^s aa [Hj. Thus, one expects that the production 
pattern of heavy quarks is the same as light ones, once the collision energy is high 
enough. 

Unfortunately, there is no calculation of the heavy quark production cross section 
for the AA reaction. Therefore, to get a numerical results which can be compared 
with RHIC experiments one uses phenomenological models which take into account 
the important features of the heavy quark production by strong color fields. This 
can be done by expanding exponents in Eq. (J2J) (together with terms describing gluon 
rescattering omitted here) to the leading order in Q 2 and Fourier-transforming to the 
momentum space. This way one gets the cross section for the quark production in pp 
collisions in fc^-factorized form [7] 

d 2 pdy 1 dy 2 = / ^ / ^ 92 ^ Xl '^ *' "'ii'^) <^ X2 ' £3)' ( 5 ) 

where A is the amplitude for the subprocess process g*(q 1 )g*{q 2 ) — * 9.999- It reduces 
to the familiar perturbative amplitude for on-shell gluons gg — > qq after taking the 
limit 9 2 ,9 2 and averaging over azimuthal angles of q 1 and q^. However, this 
approximation turns out to be not valid at high energies since it disagrees with the 
experimental data on charm and bottom production at Fermilab. 4>(x, q) in (SJ is 
the unintegrated gluon distributions such that <j)(x,q) = dxG(x,q)/dq 2 . In principle 
<j>(x,q) can be found from the QCD evolution equations at high energy 0. However, 
finding the numerical solution to those equations is still not completely solved problem 
even in the mean-field approximation. We rather prefer to use a simple model for the 
unintegrated gluon distributions based on available theoretical information about the 
solution to the evolution equations. Assume that the cross section for the quark 
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Figure 1. Dependence of the charmed meson yield on centrality at central and 
forward rapidities for dAu and AuAu collisions. 



production in pA and AA reactions can be written in the same kx factorized form as 
for pp case (JSJ). Then, (p(x,q) must satisfy the following conditions: 

a s 2 (2tt) 2 



( x >i) = Z o^,4 HQl/f) , i < Ql saturation (6) 



0(x,g 2 )~g 2 , q 2 > k 2 gcom , pQCD behavior (7) 
where fcg Com 3> Q 2 is the "extended geometric scaling" [§] momentum squared, and 
<p(x,q 2 ) = <p(q 2 /Q 2 s ) , q 2 < fc 2 com , "extended geometric scaling" (8) 

Of course there is more than one function satisfying these conditions. In my paper 
with D. Kharzeev [H] a parametrization is suggested which is inspired by the double 
logarithmic approximation to the DGLAP and BFKL equations. Our numerical 
calculations have been done with that function. The results of calculations are 
presented in Fig. ^ an d Fig. 

The charmed meson yield shown in Fig. ^ gets suppressed as one goes from 
central rapidity y = to forward rapidity r\ = 2 both in pA and AA collisions. 
This is because the saturation scale of a nucleus increases with rapidity. At RHIC 
Q 2 ~ 2e 3y GeV 2 . Therefore, at midrapidity Q 2 < m 2 and the charmed quarks are 
produced perturbatively. As the result yield scales like 
dN AA , 
dn '*)=° 

At forward rapidity Ql > m 2 implying that 

dNAA, 
dn 

The same formulas hold for dA after replacement of 7V part by N^ t and A co n by -/V^jj. 

The suppression effect is insensitive to high energy quantum evolution effects 
which change the anomalous dimension of the gluon distribution functions for 
transverse momenta satisfying Q 2 < q 2 < fc 2 com . If however one applies an infrared 
cutoff on a spectrum at transverse momenta of the order of Q s , then the region 
q 2 < fcg Com where the quantum effects are important will dominate the yield. The 

resulting nuclear modification factors scale like Raa ~ I /A 1 / 3 ~ l/N^ t ~ 0.5 and 



A 4 / 3 ~ iVcon . (9) 
implying that 

A - AW ■ (10) 
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Figure 2. \(D° + D°) spectra in heavy-ion collisions at i/s = 130 GeV and 
y/s = 200 GeV. 



RdA ~ l/A 1 / 6 ~ l/JN^ ~ 0.75 in the forward region. As we have already noted, 

these suppression factors at forward rapidities are of the same order of magnitude 
as in inclusive meson production case, since at low- re (forward rapidities) production 
pattern of charmed quarks is the same as the light ones. 

In Fig. |21 we compare the open charm spectrum obtained from our saturation 
model based on fey-factorization with predictions based on collinear factorization 
approach (PYTHIA). We observe that PYTHIA spectrum is significantly softer than 
ours. Therefore, we predict much harder open charm spectrum. Of course, the 
steepness of the PYTHIA spectrum depends on the intrinsic momentum parameter kg . 
Our spectrum can probably be reproduced by PYTHIA if one takes ko ~ Q s » Aqcd 
and fixes the i^-factor. This large value of intrinsic kr would however signal the 
breakdown of the collinear factorization approach. 
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